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ABSTRACT 

We consider the effect of a long range, flavor changing tensor interaction 
of possible gravitational origin. Neutrino mixing experiments provide the 
most sensitive probe to date for such forces — testing the equivalence prin- 
ciple at levels below 10~ 20 . Here we justify and generalize a formalism for 
describing such effects. The constraints from neutrino mixing experiments on 
gravitationally induced mixing are calculated. Our detailed analysis of the 
atmospheric neutrino data confirms a remarkable result: the atmospheric 
neutrino data implies the same size force as does the solar neutrino data. 
Additional tests of this suggestive result are discussed. 
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1 Introduction and Overview 



All experiments sensitive to solar neutrinos have measured a flux of elec- 
tron neutrinos smaller than that predicted by various solar models]]]], [[J, 
[f|]. It is widely accepted that the solar neutrino deficit is a result of neu- 
trino flavor mixing. The popular mechanism for this mixing is ascribed to 
nondegenerate neutrino masses and the inequivalence of the neutrino weak 
and mass eigenstates^J. A great many studies have been devoted to deriv- 
ing constraints on the mixing parameters: Am 2 , the neutrino mass-squared 
difference, and sin 2 29, the mixing angle, from the solar neutrino data. 

In this paper we consider an alternative mechanism of neutrino mixing in 
which the flavor changing interaction is assumed to be gravity ||, ||. This 
assumed violation of the principle of equivalence by neutrinos will lead to 
flavor mixing, even if the neutrinos are massless. Consequently, data from 
neutrino mixing experiments can be used to test how well the equivalence 
principle is obeyed by neutrinos. As will be seen, current data probe the 
equivalence principle to a suprisingly stringent level. 

In addition to solar neutrino experiments, recent studies of atmospheric 



neutrinosO, M, M, 1 10 ] also point to the existence of neutrino mixing. It is 



rather intriguing to note that, in contrast to the conventional mass mixing 
mechanism, the gravitationally induced mixing mechanism can account for 
the combined solar and atmospheric neutrino data by assuming only v e — 
mixing. In the remainder of this paper we shall discuss the phenomenologi- 
cal consequences of such a flavor dependent gravitational interaction of the 
neutrinos. 
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We begin in Section II with a discussion of how neutrinos probe for flavor 
dependence in long range forces. We start with the time delay in a gravita- 
tional field for a massless particle and its previous applications to photons 
and neutrinos. Next is a discussion of how this same time delay, if flavor de- 
pendent, can produce neutrino oscillations. Then our formalim is presented; 
first for massless neutrinos and then for massive neutrinos. When nonzero 
neutrino masses are included into our formalism, we find a new, physically 
significant phase parameter. We conclude Section II with a discussion of the 
ambiguities of the formalism, and how some of them can be avoided. In 
Section III we summarize the constraints from accelerator and solar neutrino 
data; and we present a new result for the atmospheric neutrino anomaly. 
Also, the parameter regions that planned experiments can probe are esti- 
mated. In Section IV we compare the allowed parameter regions obtained in 
Section III and find the remarkable result that atmospheric and solar obser- 
vations both select the same, small parameter region for equivalence principle 
breaking. We discuss more stringent tests of this radical idea and mention 
some possible outstanding theoretical questions that it generates. 

For those readers disinclined to move beyond the orthodoxy of metric 
theories of gravitation, which preserve the equivalence principle, the results 
presented here can be used to set the most restricive bound to date on the 
validity of the equivalence principle. They can also be reinterpreted as ev- 
idence for the existence of a new, very long range, flavor changing tensor 
interaction of gravitational strength. 
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2 Foundations and Formalisms 



2.1 Time Delay of Massless Particles in a Gravita- 
tional Field 

The most familiar tests of the weak equivalence principle (WEP) are exper- 
iments of the Eotvos-type ]l 1|| which measure the gravitational acceleration 
of macroscopic objects. Currently, it is found[Q that gravity accelerates all 
macroscopic objects at the same rate to an accuracy of one part in 10 12 . 

A test of the WEP for massless particles stems from Shapiro's obser- 
vation that the transit time for a photon traversing a given distance in a 
gravitational potential, 0(r(t)), is dilated by the amount ]TB 



At = -(l+ 7 ) / <Kr(t))dt, (1) 

where 7 is a parameter in the parametrized post-Newtonian (PPN) formalism 
[|T4]j and the integration is from the time of emission of the photon to when 
it is detected. In general relativity, 7 = 1 and radar ranging experiments [|l5l 
have verified this prediction. In alternate theories of gravity encompassed by 
the PPN formalism and satisfying the WEP, 7 need not equal unity, but it 
must be the same for all particles. Thus, observation of different 7-values 
for two different particle species would constitute a violation of the WEP. 
Limits on such a violation have been obtained from supernova SN1987A by 
comparing the arrival time of photons and neutrinos with the result [|16|1 , |J 



|7 7 -7^| < few x 1(T 3 . (2) 
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In the same way, comparing neutrinos with antineutrinos yields ||IB|| [IS 

|7, e -7pJ<10- 6 . 



(3) 



2.2 Neutrino Oscillations from Flavor Dependent 7- 
Values 

As a prelude, we remind the reader of the famous COW experiment p0| . In 



this experiment, a beam of neutrons is passed through a beam splitter and the 
separated beams are sent through two spatial regions that are gravitationally 
inequivalent. Through this passage the two beams acquire a phase difference 
determined by integrals of the gravitational potential through which they 
pass. The split beams are reconstituted and the diffraction pattern resulting 
from the gravitationally acquired phase difference is observed. 

Similary, if 7^ is flavor dependent, then different neutrinos will undergo 
different gravitational time delays when passing through the same gravita- 
tional potential and thereby acquire different phase shifts. These phase shifts 
are observable owing to the difference in the particle bases that diagonalize 
the weak and the gravitational interactions. For instance, in the case of two 
neutrino flavors, the weak basis (u e , v^) may be related to the gravitational 
basis (ui, v-i) as 



(4) 



where Qq is the mixing angle. As a consequence, a v e will be able to oscillate 
into v^. (It should be stressed that this oscillation will take place even if the 
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neutrinos are massless.) More specifically, assuming plane wave propagation 
and a constant potential, the phase difference due to the time delay acquired 
by neutrinos of energy E traversing a distance / is 

S t = (7i - -ftWE. (5) 

This phase shift was first described by Gasperini[|| in terms of a nonuni- 
versal gravitational red shift. However, this result does not include the contri- 
bution of phases from the spatial part of the wavefunction. In order to com- 
bine these two effect in a systematic way, a more specific model is required. 
Such a model was considered for spinless neutrinos]^]. For completeness, we 
extend that discussion explicitly to the spin-1/2 case belowfl2l|. 

In the absence of non-gravitational interactions, the properties of a spin- 
1/2 particle in a specified gravitational field, G a p, are usually described to 
first order (linearized therory) by the interaction Lagrangian density 122 



Ant = i f -G a ^ la d^ - {dcfrnf&X, (6) 

where / = ^JS-kGn, Gat is Newton's constant and the metric of flat space is 
g a p = (+1, —1, —1, —1). For our purpose, we postulate an interaction of the 
above form but one which allows the neutrinos v\ and v% to couple to gravity 
with different strengths J\ and ji- 

This breakdown in universality of the gravitational coupling strength de- 
stroys the symmetry that keeps the graviton massless Currently, no 



attractive theories have been proposed which can break the equivalence prin- 
ciple and yet keep the graviton massless in order to reproduce the experi- 



mental result for the deflection of light 1 24]. However, such a theory may 



be possible |2lJ, and it has been shown that a theory with mixed massless 



and massive gravitons can be made consistent with all the observations [p€ 
Since the purpose of the model discussed here is only to put together the neu- 
trino's temporal and spatial phase shift differences, and since the graviton 
mass is not directly relevant to the neutrino phenomenology, we put aside 
the formidable question of the complete consistency of such a theory. 

The postulated interaction leads to the equations of motion for the mass- 
less neutrino fields, Uj, 

+ f -lG^) la d p + k{a a( y#) lfl \ Vj = 0, j = 1, 2, ... (7) 

We limit our discussion to situations in which G a ^ varies very slowly on 
a scale of order the neutrino Compton wavelength, and therefore drop the 
terms involving a derivative of G a ^. In this case, we readily find that the Uj 
satisfy a Klein-Gordon equation, 

(g aP + fjG^dadM = 0, (8) 

where 0(G 2 ) terms have been dropped for consistency with a linearized the- 
ory. We assume the gravitational field is determined by a static macroscopic 
matter distribution in the harmonic gauge. (Since general covariance is bro- 
ken, the result will in fact be gauge dependent.) Such a field is given in terms 
of the Newtonian potential <fi by|2"7|] 

Ga.fi = 20<W/, (9) 

where d(oo) — > 0. 
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We note from Eqs. (^) and (|9|) that the presence of the gravitational field 
modifies the flat space metric in a neutrino species dependent way such that, 



for Uj, 



2f. 

9a/3 ~> (9a/3 + -rWap)- (10) 

Comparing this with Eq. (1.1) in Ref. (note that is defined to be 



positive in Ref. |28[ ) , we see that our approach amounts to the case in which 
the PPN parameters a and 7 are identical and the parameters fj are related 
to the PPN parameters by 

fj = fir (11) 

To illustrate the essential properties of the resulting phase shifts, we con- 
sider the case of constant 0, where we have the energy-momentum relation 

£ 2 (l + 2 7 ^)=p 2 (l-2 7j 0). (12) 

To first order in 0, the energy eigenvalues of the neutrinos, Uj, having the 
same momentum are given by 

E 3 = (l-2 l3 <P)p. (13) 

For the simple case of two neutrinos, this implies that, after traversing a 
distance I, the two components, (u±, 1/2), of a state v e will develop a phase 
difference of 5 = 2(71 —72)0/^- If we revert to a description of the oscillation 
phenomenon utilizing states of equal energy but differing momenta, the phase 
difference becomes, to first order in 0, 

5 = 2(7i-72)0^, (14) 
which is twice that obtained from the Shapiro effect alone. 



If we compare this phase shift with that obtained in the well known case 
of vacuum oscillations induced by a neutrino mass difference, we find that 
they are related by the formal connection, 

Am 2 



2E 



2£|0|A 7 , (15) 



where A7 = 72 — 7 i . By analogy, the u e survival probability after traversing 
a distance I is given by 

Til 

P(u e -> v e ) = 1 - sin 2 26 G sin 2 — (16) 

A 

where A is the oscillation length which is here given by 

, / 1(T 20 wlO GeVx , , 

A=6 - 2km WHH (17) 

Thus, in sharp constrast to the case of oscillations induced by a neutrino mass 
difference where A grows with energy, gravitationally induced oscillations 
are characterized by an oscillation length that diminishes with increasing 
neutrino energy. The two mechanisms may therefore be distinguished by 



measuring the neutrino energy spectrum [29], [3C], [31 



2.3 Neutrino Oscillations From Both Gravity and Mass 
Terms. 

In the preceding discussions we have described how neutrino mixing is gen- 
erated by breaking the equivalence principle. It is important to note that 
this mixing may not preclude the generation of neutrino mixing from vac- 
uum mass terms — which is the more commonly discussed source of neutrino 
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mixing. It is possible that there is flavor dependence in both the neutrino 
mass matrix and in the gravitational coupling matrix, simultaneously. This 
possibility has been discussed in the literature [Q, j3l| . However, these dis- 
cussions have overlooked an important point which we explain now. 

For simplicity, we shall assume there are only two neutrino flavors. Choos- 
ing to work in the basis which diagonalizes the charged lepton mass matrix, 
the evolution of neutrino flavor in a medium is described by 



■£ 



Am 2 



-U 



M 



-1 
1 



u 



M 



+ £|0(r)|A 7 U, 



+ ^G F N e 



G 

1 
-1 



-1 
1 



(18) 



Here Am 2 = m 2 . — m 2 denotes the difference in neutrino vacuum masses, 
N e is the electron density of the medium, and Gp is Fermi's constant. The 
first term describes the contribution from vacuum masses, the second term 
describes the contribution from equivalence principle breaking, and the third 



term describes the contribution from a background of normal matter [32 



There are two unitary matrices which parametrize the mixing, Um and Uc, 
and they are generally completely unrelated. The subscripts M and G denote 
quantities which come from the vacuum mass term and from the gravitational 
mixing, respectively. A general representation for a unitary matrix is given 
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by 



U = e lx 



-iff 







(19) 



: l0i cos# sin i 

e ia - sin 9 cos 9 J 
where Xi a an d P denote arbitrary phases and 9 denotes an arbitrary an- 
gle. However, not all of these phases are observable in neutrino oscillation 
experiments. 

As is well known, when there is no gravitational mixing (e.g., in a vanish- 
ing gravitational field), all the phases in U M may be eliminated by redefini- 
tion of the spinors, e.g., v' e = e l2aM u e . Then only 9 M is observable. Similarly, 
when the contribution from the neutrino vacuum masses is negligble (e.g., 
in a large gravitational field), all of the phases in U G may be eliminated and 
only 9 G is observable. Thus 9m and 9 G are each independently observable 
parameters. Notice that in general, when neutrino mixing receives contribu- 
tions from both the gravitational and vacuum mass terms, the x' s an d /9's 
can be eliminated but the a's cannot. Putting all of the residual phase into 
the gravitational term, the flavor evolution may be parametrized as: 

d 
i— 
dt 







Am 2 


— cos2^ M sin26 ) M 


{ "e ) 




AE 


sin 29m cos 29 m 



+ £|0(r)|A 7 



+ ^G F N e 



- cos 29 G e~ t2a sin 29 G 
e i2a zm29 G cos2# G 



(20) 



1 












1 




-1 


> 







where a = (a G — ckm)- Since this phase cannot be eliminated by redefinition 
of the spinors, it will have experimental consequences. 
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In the case of constant backgrounds (iV e and 0), this equation can be 
easily solved, yielding the oscillation probability 

P(u e , x - Vfl , y) = I j^TM { X ~ cos ( 2 ^+^(y ~ *))} > ( 21 ) 

where b is the diagonal element of the total mixing matrix 

Am 2 \/2 
b = — — cos20 A / + £|0|A7cos20 G - —G F N e (22) 
AE 2 



a is the magnitude of the off-diagonal element 

Am 2 



■ sin 20 M + E\<f>\ A7 sin 26» G e 



-i2a 



(23) 



4£ 

and x and 7/ denote the position of the neutrino source and detector, respec- 
tively. This solution explicitly shows that the phase is only observable when 
both neutrino vacuum mass mixing and gravitationally induced mixing are 
relevant, i.e., when Am 2 , A7, Qq and 9m are all nonzero. 

Phases in the vacuum mass mixing matrix may lead to violations of time 
reversal (T) symmetry. Under this symmetry, neutrino oscillation probabili- 
ties transform as|33| 



P(vp,x -+v e ,y) — ► P(u e ,y -> v»x). (24) 

Unitarity for two neutrino flavors implies that 

P(i/ M , x -> u e , y) = P(y e , x -> i/p, y) (25) 

It follows that a difference in the probabilities P(v e , x — > u^, y) and P(u e , y — > 
z/^, x) is a measure of T violation. However, note that the presence of a back- 
ground medium may lead to an apparent T- violation, even though the under- 
lying dynamics is T-invariant. To avoid this one must consider a background 
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medium which is symmetric about the midpoint of neutrino propagation, 
e.g., a constant background as in Eq. (pll). It is then clear that the phase a 
does not lead to observable T-violations in two-neutrino oscillations. 

In summary, for two neutrinos there are two quantities which parametrize 
the strength at which mixing occurs, A7 and Am 2 , and three quantities 
which parametrize the mixing, 9q, 9m and a. Previous analyses omitted the 
phase a, but in a general analysis it should be included. We plan to study 
the effects of the phase a in more detail in a future communication. 



2.4 Potential Ambiguities. 

It is necessary to choose a metric in order to confront experiments. For 
any metric choice, breaking the equivalence principle implies neutrino mix- 
ing. However, because the particular breaking of the equivalence principle 
considered here also breaks general coordinate invariance, physical results 
depend on the choice of metric used. At present, there are few hints as to 
Nature's choice. 

In deriving our formalism, we assumed the harmonic metric. Although 
this is commonly done in most discussions of K—K oscillations [Q, where one 
can test the equivalence principle between a particle and its antiparticle, the 
harmonic metric is certainly not the only possible assumption. For example, 
the authors in Ref. |21| assumed the Schwarzschild metric instead. They 
have obtained a neutrino oscillation formalism which only differs from ours 
by parameter redefinitions. 

Some metric choices do, however, lead to observable consequences. For 
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instance, the oscillation wavelength in Eq. ( |TT| ) is not invariant under shifts 
of by an overall constant. Thus the expression for that should be used 
in calculations is ambiguous, but this ambiguity must be resolved in order to 
derive phenomenological constraints. The most common procedure for fixing 
a definite value for is to assume that it vanishes at an infinite distance 
from the source. This insures that as gravity is "turned off" the results of 
special relativity are recovered. With this assumption, the local gravitational 
potentials due to various sources are given in Table 1. We see that the value 
of this potential anywhere in our solar system is dominated by the mass 
distribution on scales larger than the galaxy. However, since the details of the 
structure on extra-galactic scales is still a subject of investigation, the precise 
value of from these scales is somewhat uncertain. But it is apparent that, 
because the dominant comes from scales much larger than the propagation 
length in any current neutrino observation, it is reasonable to ignore the 
variation of over the solar system and to take it to be a constant — the 
same constant for all local neutrino experiments. Since A7 and always 
occur multiplied together, we can avoid the issue of the uncertainty in by 
calculating the constraints on the product 0A7. Calculating the constraints 
in this manner has several additional advantages. For the most part, it covers 
the possibility that does not "turn-off" at large distances, or that there is 
some " cosmological" value for 0. Also, it covers the possibility that is 
generated by a long range tensor interaction other than gravity. Thus, by 
treating as an unknown constant for our solar system, we avoid many of 
the ambiguities inherent to these calculations. 

We briefly mention the interesting possibility that the coupling between 
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gravity and neutrinos might be spatially anisotropic. Such effects have been 
considered for photons (see, e.g., |35j and references therein). For example, 
the equivalence principle breaking might come from a coupling of neutrinos 
to V0, the gradient of the gravitational potential, instead of to <fi. This 
coupling exists in our formalism but its effects are expected to be drastically 
suppressed compared to the leading effects considered here (see the discussion 
after Eq. (|7])). This kind of coupling may also occur in string theories^ 



Constraints on such a coupling have been considered in Ref. |30| . An inter- 
esting consequence of this coupling is that the amount of neutrino mixing 
depends on the orientation of the neutrino's momentum with respect to the 
gravitational source. 



3 Experimental Constraints. 

In this section we calculate the experimental limits on the gravitationally 
induced mixing (gim) parameters, A7 and 8g, from neutrino experiments. 

3.1 Constraints from Rare Muon Decays. 

The anomalous neutrino-gravity interaction breaks the lepton number sym- 
metry of each family and will give rise to rare leptonic decay processes. 
One might expect a particularly stringent constraint to arise from the de- 
cay /i — > cy. Here we consider what limits are placed on the gim parameters 
by this process. 

The specific gim interaction considered in the preceding section affects 
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only neutrinos. Consequently the leading contributions to the decay /i — > e 7 
arise at the one-loop level. These contributions turn out not to be renormaliz- 
able. However, it might be argued that the equivalence principle is broken at 
a scale that still preserves the weak isospin symmetry. Since gravitation is not 
(yet) describable as a renormalizable gauge theory, the argument is suspect. 
Nevertheless, we shall examine this scenario, since it will generate the rare 
muon decay at tree level and presumably provide the most restrictive calcu- 
lable constraint on the equivalence principle violating parameters obtainable 
from this process. While the Standard Model symmetry allows an additional 
equivalence principle violating term for the right-handed charged leptons, we 
shall ignore this complication. We simply replace the (left-handed) neutrino 
fields in their anomalous gravitation interactions by the corresponding left- 
handed family isospin doublets and make the derivatives appearing there 
SU(2) gauge covariant in the usual fashion. This leads to three tree-level 
Feynman diagrams which yield a decay rate of order magnitude 

T ~ am M (A70 sin 26 G ) 2 , (26) 

where a is the fine structure constant and m M is the muon mass. Taking 
the large but somewhat uncertain value of \<ft\ = 3 x 10~ 5 , we estimate the 
branching ratio for fi — > e 7 to be 

B(n -> erf) ~ 2.3 x 10 6 (A 7 sin 29 G ) 2 , (27) 



Since the current experimental limit[37 is B(/j — > e 7 ) < 5 x 10 n , this leads 
to the constraint 

|A 7 sin2fl G | < 5 x 1(T 9 . (28) 
16 



Except for very small mixing angle, this constraint on A7 is inferior to the 
constraints obtained below from solar, atmospheric and accelarator neutrino 
experiments. 

3.2 Constraints from Accelerator Experiments. 

There exist numerous constraints on neutrino mixing from experiments using 
neutrinos produced by accelerators. Because accelerators produce a neutrino 
flux which is mostly muon-neutrinos with a much smaller electron-neutrino 
component, these experiments mostly probe — u x mixing with a reduced 
sensitivity to u e — u x mixing. Accelerator experiments typically involve neutri- 
nos with energies ranging from a GeV to hundreds of GeV, and propagation 
lengths up to a kilometer. Since the neutrino cross section increases with 
energy, these experiments can achieve large event rates and so can probe for 
relatively small values of the mixing angle. Most importantly for us, the high 
neutrino energies in these experiments make them powerful tools for testing 
the equivalence principle because the gravitationally induced mixing 'turns 
on' with increasing energy. 

It is difficult to utilize the published analyses to obtain the exclusion 
region for I0A7I and sin 2 20 G . This is because the energy dependence of gim 
is so radically different from the case of mass mixing that the analyses must 
be completely redone. However, many important details such as the energy 
dependence of the observed and expected event rates are seldom published. 
We are thus forced to a less than optimal estimate of the exclusion region. 
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Here we calculate the exclusion region for acclerator experiments in the long 
wavelength and the short wavelength limits. These lines are then extended 
to intersect in order to estimate the full exclusion region. The constraint 
in the short wavelength limit is obtained by taking the average probability 
at minimum sensitivity to be the same for both mixing mechanisms. The 
constraint in the long wavelength limit is obtained by using the minimum 
Am 2 from the published analysis to calculate the typical energy scale of the 
experiment and, with the average length given in the published reference, we 
compute the minimum value for |0A7|. In this admittedly crude fashion, we 
obtain the excluded regions shown by the straight-line contours in Fig. (1). 
We estimate that these contours are accurate up to factors of 3 in |0A7|. 
We have selected experiments ||38|| having the largest values of E ■ I, for they 
provide the most stringent limits. 

Recently, the LSND accelerator neutrino experiment has reported some 
candidate events that may||39||, or may not [[5]], indicate — v e neutrino 
oscillations. The favored parameter region lies above |0A7| > 5x 10~ 18 and so 
large mixing angles are excluded by other, higher energy, accelerator neutrino 



experiments |[41||. There may be a very small parameter region, allowed by all 
accelerator experimental results, at sin 2 29 w 10~ 3 and |</»A7| ~ 2 x 10~ 16 . 
In light of the preliminary nature of the LSND experimental results, we shall 
not pursue this analysis further. 

There are several serious proposals for a new generation of accelerator ex- 
periments with neutrino propagation lengths of several hundreds kilometers [[42] 
43[| . For example, the MINOS experiment would send a neutrino beam from 



the Fermilab accelerator to the Soudan-2 detector located 730 kilometers 
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away. The lower, curved contours in Fig. (1) estimate the gim parameter 
region that can be probed by this experiment. Note that matter effects start 
to become important at these longer distances and cause a difference in the 
neutrino and antineutrino oscillation probabilities, as shown by the solid and 
dotted curves in Fig. (1). 

3.3 Constraints from Reactor Experiments. 

Numerous constraints on neutrino mixing have been obtained from experi- 
ments using neutrinos produced by nuclear reactors. Nuclear reactions cre- 
ate neutrinos with energies ranging from a fraction of an MeV to order 10 
MeV, and these have been detected hundreds of meters from the reactors. 
While these relatively low energy neutrinos provide stringent constraints on 
mass mixing, they provide only very weak constraints on gim because of the 
contrasting energy dependence of the two mechanisms. Specifically, reactor 
experiments only probe values of |0A7| > 1 x 10~ 16 . For the gim mecha- 
nism, reactor constraints are completely surpassed by those from accelerator 
experiments, so we shall not discuss them further. 



3.4 Solar Neutrino Constraints. 

The nuclear reactions which power the sun also produce neutrinos. These 



neutrinos have been observed here on Earth by several experiments! 44] [45 



[ |46|j |f47| . These experiments study different ranges of the neutrino spectrum, 
which extends up to 14 MeV and is a superposition of several components 
(see, e.g., |l|). All of these measurements find far fewer neutrinos than ex- 
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pected (see Table 2). This is commonly interpreted as evidence for neutrino 
mixing. This is because both the production and detection of solar neutrinos 
primarily involve only the v e flavor, so any mixing will reduce the observed 
flux. Here we shall assume that all of the neutrino mixing comes from the 
gim mechanism, and derive the implications of solar measurements for gim 
parameters. 

The oscillation probabilities are calculated analytically. Here, in accor- 
dance with the discussion in Section 2.4, we shall take to be a constant (in 
our earlier analysis[|29| we used the generated by the Sun). Gim effects start 
being important for solar neutrinos when |0A7| > 2 x 10 -25 . This is when a 
10 MeV neutrino undergoes half of an oscillation in its propagation from the 
Sun to the Earth. For larger values of |</>A7| the oscillation wavelength, Eq. 
(|TTD, is smaller, and when the wavelength is less than the size of the Sun, the 
effects of the background matter become important. The analysis of matter 
effects [32| parallels the well known case for mass mixing[Q (for a review, 



see, e.g., |4"5||), with the substitutions of Eq. flT5| ) (for constant potential) 
and 9 — > 6q. The condition for a resonance to occur is now given by 

V2G F N e = 2E\(f>\ A 7 cos 29 G , (29) 

where N e is the electron density in the sun. Note that when matter effects 
are relevant, the sign of A7 is important. For A 7 > the neutrinos can go 
through a resonance, and if the transition is adiabatic, large reductions in 
the neutrino flux occur. The calculated oscillation probability, for A7 > 0, is 
shown in Fig. (2) as a function of neutrino energy, E, times \<p\ A7. Scanning 
across Fig. (2) from low to high energies, we see that first the long wavelength 
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oscillations occur, then the nonadiabatic side of the MSW well, and then 
the adiabatic side. This order is a reversal from that of the mass mixing 
mechanism, which provides a means to experimentally distinguish the two 
mechanisms. 

The most recent solar neutrino data and solar model results are used in 
our calculation of the constraints on gim parameters. In addition to the re- 
sults given in Table 2, the Kamiokande group has looked at how their data 



depends on energy ||45|| . These data can be used to constrain mixing, since 
mixing effects can be energy dependent (see Fig. (2)). Much of the neutrino 
energy dependence is lost in the neutrino-electron scattering process, which 
occurs in the Kamiokande detector, but enough remains to be important. We 
have folded our oscillation probabilities into the neutrino spectrum and then 
integrated this times the cross section over energy to get the predicted result 
for each Kamiokande energy bin, and for each of the other experiments. The 
total chi-squared is calculated between the predicted and measured results 
and this is used to find the allowed gim parameter regions. Only the experi- 
mental errors were included, since these are estimated to be the largest. We 
first calculated the allowed parameter region in the mass mixing mechanism. 



We then compared these calculations with those of others (e.g., |5(|) and 
found that our 90% contours were in good agreement with theirs. We used 
the same value of chi-squared to calculate the corresponding 90% contours 
in the gim mechanism. This insures that the mass mixing mechanism and 
the gim mechanism are treated equally in fitting the same data. 

The results of our analysis are given in Fig. (3). There are plotted 
the allowed parameter regions at 90% and 99% confidence levels, assuming 
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A7 > 0. For A7 < there is no allowed parameter region because then 
matter effects suppress mixing. The updated data, and the use of a constant 
|0|, do not significantly alter the allowed parameter regions from those found 
in our earlier analysis |29[, which were independently confirmed[[3(J. There are 
two allowed parameter regions: one at large mixings, where the MSW effect 
suppresses mostly the intermediate and high energy neutrinos, and another 
at small mixings, where the MSW effect suppresses mostly the intermediate 
energy neutrinos. These two parameter regions are analogous to the two 
MSW regions found for the mass mixing mechanism (see, e.g., [jjOfl ). 

For the mass mixing mechanism there is also a long wavelength, vacuum 
oscillation solution (see, e.g., |5l|). For the gim mechanism, long wavelenth 
effects are relevant for large vacuum mixings and I0A7I ~ 2 x 10~ 25 . Around 
this point only the highest energy solar neutrinos are suppressed (see Fig. 
(2)). However, the data require some suppression at low neutrino energies so 
there is no allowed gim parameter region near this value [p9fl . 

There are several new solar neutrino experiments which will test these 



explanations of the data. The SNO[|52| experiment will definitively test all 



neutrino mixing explanations of the solar neutrino deficit by performing a 
flavor independent measurement of the solar neutrino flux. The SNO and 
Super-Kamiokande||53|| experiments will be able to measure the energy de- 



pendence of the high energy solar neutrino flux. The large mixing solution 
would not give an energy dependence observable in these experiments, but 
the small mixing solution would. Both of these experiments can also look for 
day-night variations from a resonance with matter in the Earth. These varia- 
tions are large for these experiments when the mixing angle is large and when 
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0.75 < |0|A7/1(T 20 < 2.5. Thus, if the experiments are sensitve to small 
day-night variations, part of the large mixing solution may be probed by this 
method. In addition, the BOREXINO experiment ||54|| will measure the flux 
of Be 7 neutrinos with high statistics. This flux is almost monoenergetic so 
long-wavelength oscillations are important for a wide range of parameters. 
By looking for temporal variations on scales from a few days to a few months, 
correlated with changes in the Earth-Sun distance, this experiment can probe 
the parameter region 



3 x 1CT 4 < 1^ < 0.3 (30) 
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at large values of the mixing angles |p5|. Thus BOREXINO can probe all of 
the large mixing angle solution. These new probes using neutral current flux 
measurements, spectral distortions and temporal variations are especially 
important because they are independent of uncertainties in the solar model. 



3.5 Connection with Supernova Dynamics. 

In a stellar collapse, neutrinos are created in large numbers in the hot, high 
density core, and then diffuse out to lower matter densities where they then 
freely pass through the outer layers of the star. If the neutrinos were to go 
through a resonant flavor conversion in the region near the supernova's core, 
this would have important, observable effects. Present models of supernovae 
suggest that this would increase the size of the explosion by about 50%, and 
it would also block the production of r-process nuclides [|5(J. The first effect 
would be a welcome one, since historically theoretical models have fallen 
short of the observed energy. However, the second effect would eliminate 
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one of the most promising sites for r-process nucleosynthesis (see, e.g., [p7[ 
and references therein). Our knowledge of supernova dynamics is presently 
too sketchy to allow reliable constraints to be placed on any neutrino mixing 
parameters. Here we estimate the gim parameter region that is relevant for 
supernova dynamics. 

Qualitatively, supernova neutrinos propagate from high to low densities 
just as do solar neutrinos. However, calculating the probability that a su- 
pernova neutrino's flavor survives a resonant transition is technically more 
difficult because the large neutrino background in a supernova makes the 
flavor evolution nonlinear [5S|. A recent numerical study has found that this 



effect substantially reduces the relevant parameter region|59[ . Using the re- 



sults of this study, we estimate that gim parameters in the range 

A 7 >4xl0" 14 , sin 2 2#>l(r 3 (31) 

would have substantial effects on supernova dynamics. In this estimate we 
took the gravitational potential at the nucleosynthesis region of a supernova 
to be lO" 1 . Because this value is larger than those in Table 1, we did not 
take to be a constant determined by external sources in deriving Eq. (j3l"l) 
(see the discussion in Sec. 2.4). To compare the parameter region relevant 
for supernovae with those obtained in local neutrino observations, we must 
multiply the above limit on A7 by the value of <fi relevant to the solar system. 
Taking this to be 3 x 10~ 5 , the supernova parameter region lies just above 
the parameter region probed by the solar neutrino measurements. 
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3.6 Atmospheric Neutrino Constraints. 



Recent experimental studies M M M [1C] of neutrinos produced in the Earth's 



atmosphere have found evidence for neutrino mixing. Atmospheric neutrinos 
have properties very different from solar neutrinos. Atmospheric neutrinos 
are more energetic, with energies ranging from fractions of a GeV to several 
thousand GeV. Furthermore, their propagation length is shorter; it varies 
from 20 to 13 000 km. In the mass mixing mechanism, the oscillation phase 
depends on length/ energy, so this quantity is much smaller for atmospheric 
neutrinos than it is for solar neutrinos. Consequently the mass mixing mech- 
anism requires mixings of different neutrino pairs to fit the combined solar 
and atmospheric neutrino data (see, e.g., f60fl). In contrast the oscillation 
phase in the gim mechanism depends on length ■ energy, which has the same 



order of magnitude for both the solar and the atmospheric neutrinos 2S 



This offers the suggestive possibility that the solar and atmospheric neutrino 
data can be simultaneously explained with the same neutrino mixing param- 
eters. This expectation is indeed confirmed by the results of our analysis 
discussed below. 

The most precise measurements of the atmospheric neutrino flux have 
been made with the large water Cherenkov detectors KamiokandeJ7[ |8| and 
1MB ||. These detectors each have exposures which are more than 5 times 
larger than that of any other detector. Here we confine our analysis to their 
results. Many quantites used in the calculation of neutrino event rates have 
large uncertainties. For example, the prediction for the absolute atmospheric 
neutrino event rates has an estimated uncertainty of about 30%. Thus, one 
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must be careful about which experimental results are used to constrain neu- 



trino mixing. We do not use the experimental results of FREJUS|6T| because 
of their smaller statistics, and because their results and those in Ref. |62 
depend on the neutrino-rock cross section for which sizeable corrections have 
been calculated recently^]. The predictions for the relative rate of events 
starting in the detector is believed to be much more reliable, with an es- 
timated uncertainty of only about 5%. Neutrino mixing could change the 
relative number of to v e in the flux from the expected value of approx- 
imately 2. Thus, the relative rate can be used as an indicator of neutrino 
mixing effects. 

Kamiokande and 1MB have both measured the ratio of to u e for events 
fully contained in their detector. Dividing this number by the ratio predicted 
by a Monte Carlo analysis gives 

~ Nif c /NM C ' 1 ' 

In the "sub-Gev" energy range (0.3 - 1.3 GeV for Kamiokande, and 0.3 - 1.5 
GeV for 1MB), this ratio is found to be 

R = 0.60 ±0.05 ±0.05 (Kam) 

R = 0.54 ±0.05 ±0.12 (1MB) (33) 

Kamiokande has also measured this ratio for higher energy (1.3 GeV to 
around 10 GeV) neutrinos, which involve fully contained as well as partially 
contained events. They find this "multi-Gev" ratio to be 

R = 0.57 ±0.07 ±0.07 (Kam) (34) 
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These consistent departures from 1.0 suggest neutrino mixing. 

In addition, the Kamiokande group has studied how their data depends 
on the zenith angle. As the angle between the neutrino and the zenith 
varies from to tt, the neutrino propagation length varies from 20 to 13 000 
km. Both mass mixing and gim predict that mixing effects should appear 
at "long" distances, so an unusual angular dependence could indicate neu- 
trino mixing. The cross section in the "sub-Gev" energy range is relatively 
isotropic, so no evidence for an angular dependence is expected or found 
in those data. In the "multi-GeV" energy region, the cross section is more 
directional, with the average angular spread between the neutrino and its 
associated charged lepton being rms ~ 15° — 20°. There Kamiokande has 
found an angular dependence in their data, with longer path lengths associ- 
ated with smaller R values ||. This also suggests neutrino mixing. 

Neutrino mixing modifies the observed event rates as given by 

N e = N™ c [P{v e ->v e ) + rP{u^u e )] 

N, = Nf c [P{v e ^v,)/r + P{v^v,)] (35) 

Here N a denotes the number of charged leptons of type a observed in a 
particular bin, and the MC superscript denotes the number expected without 
oscillations from a Monte Carlo calculation. The ratio of the flux to the u e 



flux is denoted by r and is approximately 2.1 [64[. The P's are the oscillation 
probablities averaged over the energy and length distributions relevant for a 
particular bin. In the calculations performed by the experimental groups, the 
data is broken up into an array of bins depending on energy and zenith angle. 
However these full data arrays are not published. For our calculations, all of 
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the sub-GeV data lies in just one bin, while the multi-GeV data is divided 
into 5 zenith angle bins. 

The oscillation probabilities are calculated analytically. Here, as for the 
earthbound accelerator neutrino experiments and the solar neutrino measure- 
ments, <p is taken to be a constant (see discussion in Section 2.4). For — u e 
oscillations matter effects are accounted for by using a two density (core and 
mantle) model for the earth, assuming A7 > 0, and taking 3/7 of the data 
to be due to antineutrinos. The oscillation probabilities are averaged over 
neutrino energy distributions as given in references || and []65| . The average 
over length distributions is calculated as 

(P„) oc ! dnP v {6)I v (6) [ dn'P{6,(j>\0',(j>') (36) 

J J bin 

where 9 is the zenith angle, P v denotes the neutrino oscillation probabil- 
ity, l v is the neutrino flux intensity at the detector without mixing, and 
P(9, 0; 8 f , cf)') denotes the probability of a neutrino with angular coordinates 
9, giving rise to a charged lepton with angular coordinates 9', </>'. For the 
sub-GeV data the last quantity is not necessary, but for the multi-GeV data 
we take it to be a gaussian with an rms spread of 20 degrees||. For I u , we 
assume that at their production point in the atmosphere the neutrino fluxes 
are independent of angle and energy. This is a good approximation for the 
multi-GeV energy range. With this assumption, dQI u at the detector is pro- 
portional to the particularly simple form of dL/L where L is the neutrino 
propagation length. 

We treat the "sub-GeV" and "multi-GeV" data sets independently, be- 
cause the relative errors (e.g., in the cross section, etc.) between these two 
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energy ranges have not been studied. For each data set we calculate the x 2 



between the predicted and observed event rates[p6[. In addition to statis- 
tical errors, we include a 30% error for the overall flux normalization. We 
do not explicitly include other, much smaller, errors in our x 2 's (e.g., those 
from particle misidentification, the Monte Carlo calculations, etc). Instead 
we first calculated the allowed parameter region in the mass mixing mech- 
anism. We then compared these calculations to those of the Kamiokande 
experimental group |§, and found that our model gave 90% contours that 
were in good agreement with theirs for particular, reasonable choices of \ 2 ■ 
We then used these values of % 2 to calculate the corresponding 90% contours 
in the gim mechanism. This procedure was adopted for several reasons: it 
tests our model and insures that it is reasonable and accurate, it includes as 
many unknown experimental effects as possible into our calculations, and it 
insures that the mass mixing mechanism and the gim mechanism are treated 
equally in fitting the same data. 

The regions of allowed gim parameters are shown in Figs. (4). Fig. (4a) 
is for v e — Ufj, mixing and Fig. (4b) is for — v T mixing with A7 > (for 
A7 < the contours are almost identical). The parameter regions allowed 
by the sub-GeV and multi-GeV data sets lie to the right of the dashed and 
solid contours, respectively. The sub-GeV measurements are insensitive to 
any angular dependence in the flux, so they allow arbitrarily large values of 
|0| A7 where the oscillations completely average out. However, the multi-GeV 
data contains some angular dependence so it extends only over a finite range 
of I0IA7. Because higher energy neutrino experiments are more sensitive to 



gim effects (see Eq. (IB) and (|ITD), the bottom of the multi-GeV contours 
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lies below the bottom of the sub-GeV contours. The overlapping region is 
where both data sets can be simultaneously explained. For v e — z/ M mixing, 
Fig. (4a), the overlapping region is relatively large and contains the best fit 
points for both data sets. Thus the sub-GeV and multi-GeV atmospheric 
neutrino data can be consistently explained by the gim mechanism. 

There are several new atmospheric neutrino measurements which will test 
this explanation of the data. Super-Kamiokande[|53[| will increase the preci- 



sion of the existing Kamiokande measurement. Soudan-2 |L0[ will measure 
the flux in a detector that is not a water Cherenkov type experiment. In addi- 
tion, the next generation of neutrino telescopes |H]], DUMAND, AMANDA, 
NESTOR and Baikal, may be able to accurately measure the atmospheric 
neutrino flux at energies which are an order of magnitude larger. While 
the mass mixing mechanism turns off at these higher energies, the gim ex- 
planation of the data does not. As these experiments measure the angular 
dependence of the flux, they will probe a parameter region an order of mag- 
nitude below that probed by the current atmospheric neutrino experiments. 



4 Discussion and Conclusions. 

We have reviewed all of the present experimental constraints on neutrino 
mixing for their implication on equivalence principle breaking parameters. 
The most stringent constraints are from: present accelerator neutrino exper- 
iments, which are sensitive down to |</>|A7 ~ 10~ 21 ; solar neutrino experi- 
ments, which are sensitive down to |0|A7 ~ 10~ 22 ; and atmospheric neutrino 
experiments, which are sensitive down to |</>|A7 ~ 10 -23 . Currently, the lat- 
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ter two types of measurements actually indicate nonzero neutrino mixing. 
Assuming that this mixing arises solely from a violation of the equivalence 
principle, we have analyzed these measurements to find the allowed values of 
the mixing parameters. A comparison of Figs. (3) and (4a) shows that the so- 
lar and atmospheric neutrino measurements can be simultaneously explained 
by v e -v^ mixing with 0.6 < sin 2 26» G < 0.9 and 2x 10~ 22 < \<f)A~f\ < 2xl0~ 21 . 
This is an extremely remarkable result since the two types of events are very 
different in their characteristic neutrino energies and propagation lengths. 

Because gravitationlly induced mixing has an energy dependence which 
is the inverse of the mass mixing mechanism, the two mechanisms give quite 
different predictions. Accelerator neutrino experiments offer a controlled, in- 
dependent test of the gim solution. These experiments already rule out the 
upper half of the parameter region allowed by the solar and atmospheric neu- 
trino data (see Fig. (1)). Next generation long-baseline experiments such as 
MINOS promise to extend these limits by over 3 orders of magnitude. This 
solution will also be tested as the next generation of solar and atmospheric 
neutrino experiments start up in the near future. In particular, the BOREX- 
INO solar neutrino experiment can use long-wavelength vacuum oscillations 
to directly probe the allowed parameter region, independent of solar model 
uncertainties. 

The test of the equivalence principle discussed in this paper is quite similar 
to the more familiar tests such as Eotvos-type experiments, Shapiro's time 
dilation experiments, and others. However, there are theoretical difficulties 
attendant with a violation of the weak equivalence principle. For instance, it 
naievely suggests the presence of a graviton mass on which there are stringent 
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experimental constraints from data on the bending of light. The importance 
of this difficulty is not clear since a consistent theory of quantum gravity has 
proved elusive. Of course, our result certainly does not compel a violation of 
the equivalence principle. An alternate interpretation might be the existence 
of a very long range tensor field, in addition to the gravitational tensor, that 
couples to electron- and muon-neutrinos in the manner described here. To 
distinguish this alternative would require positive indications with particles 
other than neutrinos — an experimentally challenging task. Neutrinos are 
uniquely suited for testing the equivalence principle because they are subject 
to only the two weakest known forces, the weak interactions and gravity. 

Experiments currently under construction will decide if the equivalence 
principle is violated, as the current evidence suggests. If this parametrization 
ultimately proves inadequate by more refined experimental results, this type 
of analysis will reinforce the concept of a geometric theory of gravity. If this 
class of mechanisms provides the proper parameterization of the data, then 
perhaps we have at hand the Balmer formula for neutrino oscillations. 
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Table 1. Values of the gravitational potential \<f>\ = G^Mjr at various po- 
sitions from various sources^]. The Sun-Sun entry is the largest value of 
1 0| in the Sun due to the Sun. The details of structure on supercluster scales 
are not well measured at present, so there is a sizeable uncertainty in the last 
entry |gj|- 



Position 


Source 


101 


Earth 


Earth 


6 x 10" 


-10 


Earth 


Sun 


1 x 10' 


-8 


Solar sys. 


Galaxy 


6 x 10- 


-7 


Solar sys. 


Virgo cluster 


1 x 10- 


-6 


Sun 


Sun 


7 x 10- 


-6 


Solar sys. 


Great Attractor 


3 x 10- 


-5 
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Table 2. Results of the Homestake^], Kamiokande-II[|3 , and combined 
SAGE [^7] and GALLEX^6f solar neutrino experiments. 



Experiment 


Process 


^threshold 


Rate (SNU) 


Theory (SNU) 


Homestake 

Kamiokande-II 

GALLEX+SAGE 


z/ e + 37 Cl^ e+ 37 Ar 
v + e — > v + e 
z/ e + 71 Ga^ e+ 71 Ge 


0.81 MeV 
7.5 MeV 
0.24 MeV 


2.5 ± 0.3 
0.50 ± 0.07* 
78 ± 10 


8.0 ± 1.0 
1.0 ± 0.15* 
137 ± 8 



*The Kamiokande flux is not given in SNU, but as a fraction of the standard 
solar model [0] prediction. 
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FIGURE CAPTION 



Fig. (1). Constraints on gim parameters from present and proposed ac- 
celerator neutrino experiments. To the right of the straight-line contours lie 
the regions of |</>A7| and sin 2 29 q that are excluded by current accelerator 



experiments (see Sec. 3.2 and Ref. [38]). To the the right of the curved con 



tours lie the regions that may be probed by MINOS, a planned long-baseline 
accelerator experiment, assuming a 10% sensitivity for a disappearance ex- 
periment and A7 > 0. The outer, solid curve is for while the inner, dotted 



curve is for v^. 



Fig. (2). Plot of P(u e -> u e ) as a function of E^A^ for sin 2 29 G = 0.4. 
The oscillation probability has been averaged over the 8 B neutrino produc- 
tion region of the Sun. 

Fig. (3). x 2 pl°t showing regions of I0IA7 and sin 2 29q allowed by the 
solar neutrino data in Table 2 at 90% (solid lines) and 99% (dashed lines) 
confidence level, assuming two neutrino v e — v x mixing. 

Fig. (4). Plot of |0| A7 versus sin 2 29q showing parameter regions allowed 
by the atmospheric neutrino data at 90% confidence level. Two flavor neu- 
trino oscillations are assumed with u e — oscillations in Fig. (4a) and — v T 
oscillations in Fig. (4b). The region allowed by the sub-GeV data lies to the 
right of the dashed contours, and the region allowed by the multi-GeV data 
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lies to the right of the solid contours. The best-fit points are also shown by 
a triangle for the sub-GeV data and a cross for the multi-GeV data. 
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This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9512220vl 



This figure "fig2-l.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9512220vl 



This figure "figl-2.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9512220vl 



This figure "fig2-2.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9512220vl 



This figure "figl-3.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9512220vl 



